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In this paper we predict many interesting new properties 
of vortices in highly anisotropic density wave systems subject 
to strong transverse electric fields. We mainly concentrate on 
ground state properties. Besides electric field-induced vortices 
we consider also thermally activated vortices. A new type of 
temperature-driven transition between two different phases 
of density waves in strong fields is predicted and several new 
properties of those phases are reported. 



I. INTRODUCTION 

It is now well realized that the dynamics of topolog- 
ical defects plays an important role in low dimensional 
electronic systems. In a particular case of density wave 
(DW) ground states (for a review see Jl|) it was found 
that topological defects are responsible for such intcrsting 
phenomena as generation of narrow-band noise Q , non- 
linear current- voltage characteristics in strong fields (the 
phenomenon associated with the so-called "phase slip- 
page" ) plH . More recently, quantum phase slip through 
creation of vortices was proposed |^] to explain low- 
temperature properties of spin-density wave compound 
(TMTSF) 2 PF 6 g§. 

This paper was originally initiated by the prediction ||] 
that DW systems should possess a novel phase in strong 
electric fields applied transverse to the direction of high- 
est conductivity. This novel phase was called "mixed 
state" in analogy with the mixed state of type-II su- 
perconductors, and was characterized by the presense of 
vortices. We have undertaken an extensive study of the 
predicted transition. The possibility of transverse electric 
field-induced transition to metallic state (type-I-like) was 
also investigated and a temperature-driven transition be- 
tween type-I and type-II regions of the phase diagram was 
predcited. We report these findings here as well as some 
interesting anticipated properties of the above-mentioned 
phases such as, e.g., non-exponential screening of electric 
fields by DW systems in " mixed state" . 



II. GINSBURG-LANDAU THEORY 

We adopt macroscopic mean-field description of the 
dynamics of DW systems and make use of (time- 
independent) Ginsburg-Landau free energy Jl| which 



is a functional of complex order parameter A(r) = 
|A(r)| exp(i0), absolute value of which is DW gap in the 
electronic spectrum. For the case of highly anisotropic 
system with an open Fermi surface: k = (±fci?(l + 
7 cos(fefcy)), k y , k z ) (here x is the direction of highest con- 
ductivity - usually the direction along conducting chains, 
y is the direction of strongest Fermi surface wrapping, 7 
is the anisotropy parameter, b is lattice constant along y; 
we neglect dispersion along z) this functional is 0: 
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We would like to emphasize the fact that the two terms 
in Finhom are of totally different origin. As is well known 
PJ, in the system described above a density modulation 
appears with the wavevector: q = qo = (2kp, ir/b, 0): 
Pspin oc Acos(qor). Spatial variation of A gives rise to 
the effective shift in q by 5q = — iV log A. Thus the gap 
is developped at the wave-vector which is slightly off its 
optimal value qo- This costs extra energy: 6E cx |<5q| - 
it is linear with <5q which is quite costly for small devia- 
tions. Instead the system rearranges its electron density 
in such a way that locally 2kp coincides with the new 
DW wave- vector. There is no more linear contribution 
to the energy, however another contribution arises - the 
energy of the non-homogenious electron density distribu- 
tion. But the fact is that this density (per chain) is: 
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and the extra energy is: SE cx p 2 oc (8q x ) 2 - it is quadratic 
in 5q x and therefore wins for small distorions of DW gap. 
This phenomenon is usually described as follows: "par- 
ticle density follows the variations of the gap". At the 
same time the gradient in y-direction does not lead to 
any change in the local charge density; its contribution 
to the free energy is of purely elastic character. 

In this paper we shall be interested in the response of 
DW system to external electric fields. We take account of 
the effect of electric field by introduction of the following 
term into the free energy functional M : 
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Here the first term describes coupling of electric field to 
the charge induced by DW distortion, the term in paren- 
thesis represents energy associated with electric field it- 
self and with the coupling of electric field to free charge 
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carriers escited over the DW-gap. The last term repre- 
sents the interation of electric field with external charges 
p[ . Strictly speaking Eqn. (0) is the action of the system 
taken with minus sign. However we shall use the term 
"free energy funciotnal" to denote a functional of which 
an extremum determines the ground state. We also em- 
phasize that electric field couples only to x-component 
of the gradient of the gap - only this component induces 
electric charge! 

Important point here is that our action is essentially 
different from that of ref. || in that it contains wave 
vector-dependent dielectric constant e = e(k). It ap- 
pears because the ground state of DW, in complete anal- 
ogy with the low-T ground state of semiconductors, is 
highly polarizable due to excitations of virtual electron- 
hole pairs. This should not be mixed neither with the 
screening due to thermally excited quasiparticles, nor 
should it be mixed with the screening due to spatial dis- 
tortion of the order parameter (see below). We stress that 
this dielectric function is generic for a gapped electronic 
system. It can be calculated using simple anisotropic 
two-band semiconductor model ffl: 
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Here Xtf is Tomas- Fermi screening length in th metallic 
state, £(fc) is anisotropic coherence length in DW state. 
(The value e(0) is usually huge: for instance e(0) ~ 10 5 
for (TMTSF) 2 PF 6 ). Now we turn to the analysis of the 
response of DW systems to external electric fields. 



III. WEAK FIELDS 

In the absense of external fields the phase of the order 
parameter A = |A|exp(i0) is the so-called degeneracy 
parameter: the free energy is globally U(l) invariant - it 
does not change if we make a global substitution 8^9 + 
const. Consequently there is a Goldstone (the so-called 
phason) gapless mode associated with this degeneracy. 
At the same time the orthogonal amplitudon mode, in 
which |A| is varied, has a gap in the spectrum at k — 
0. Therefore we expect that in the limit of weak fields 
the response of the system is given by the excitation of 
phason mode only. Thus we set | A| to be constant. After 
Fourier transform: 
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the minimization with respect to 8k becomes trivial. 
Then the remaining effective action for the electric field 
leads to an anisotropic screening length A. Screening 



length in y-direction is A y = y e{X q p)X qp and is due to in- 
teraction of electric field with quasiparticles excited over 
the gap. On the other hand screening in x-direction is 
now determined by both quasiparticles and phasons lead- 
ing to screening length X x = Xtf This is again a 
manifestaion of the fact that only the phason mode along 
x-direction couples to electric field |13| . 



IV. STRONG FIELDS 

As we have just seen, electric fields in the x-direction 
(along the chains) are screened at the short lengths 
~ Xtf at all temperatures. In contrast to that, E y com- 
ponent of the field can penetrate far enough into the bulk 
since its screening is determined be thermally excited 
quasiparticles the concentration of which is exponentially 
decreased at low temperatures. However in strong elec- 
tric fields the situation changes: A can no longer be con- 
sidered constant and it may be energetically favourable 
for the system to expell electric field from the sample by 
spatial variation of the magnitude of the gap |L4j . Firstly, 
it's obvious that in the limit E y — ► oo the state of the 
system is metallic. A simple argument for this is that 
in competition between condensate energy — n(ej?)A 2 /2 
(here n(e^) = N±/ttvf is density of states at Fermi level, 
N± is density of chains in the plane perpendicular to the 
chains) and electrostatic energy ED / 8tt the latter always 
wins in the limit of strong fields: the gap vanishes giving 
rise to screening. Thus in the region of strong fields there 
is a thin layer near the edge of the sample in which the 
state is metallic. Electric field in this layer is screened at 
the distance ~ Xtf- However, taking into account that 
actually the continuous model used here is not aplicable 
at such small distances, one should substitute Xtf for 
interchain distance. This corresponds to a metallic layer 
with the thickness of as small as a single interchain dis- 
tance. However, in the bulk the electric field is absent 
and DW state is restored. 

Thus we establish that there must exist a critical value 
of E y for which a transition occurs from DW to some 
other phase. As the field is increased from small val- 
ues two possibilities may occur: A can jump to a smooth 
configuration (Type-I) or a configuration with topological 
defects Q (Type-II). The latter means that the ground 
state is characterized by the presense of vortices - cen- 
ters the in-plane circulation of the phase of the order 
parameter around which is non-zero. These two possi- 
bilities correspond to the two types of superconductors 
in magnetic field. An ordinary isotropic superconduc- 
tor can be distinguished between these two types by the 
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evaluation of Ginsburg-Landau parameter k = A/£. Val- 
ues k < l/y2 and k > l/y/2 correspond to type-I and 
type-II respectively. In our case this simlpe argument 
is not applicable since DW materials are usually highly 
anisotropic. Indeed, usually along x-direction one has 
X x <C £ x whereas along y-direction either of X y -C £ y 
and X y S> £, y can hold depending on temperature: X y 
exponentially increases with decreasing temperature (see 
below for a more detailed discussion). Instead we shall 
use a more physical argument to decide between Type-I 
and Type-II. Namely, we shall evaluate the critical field 
for these two cases. Then we can argue that the system 
will make a transition to the state for which the critical 
field is lower. 



A. Type I 

The critical field for this case can be obtained by equat- 
ing condensate and electrostatic free energy densities: 
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Here we have introduced do — vp/2e - a parameter with 
the dimensionality of electric dipole moment (see below). 
N.B. We use the value of the dielectric constant e ~ 1 
because the metallic breakdown will occur in a thin layer 
near the boundary with the thickness ~ Xtf ~ 
polarization will not develop at such small scales! 



B. Type II - Vortices 

First of all let us clearify what is a vortex in a DW 
system. DW state is a ground state with broken U(l) 
global invariance. This means that the order parameter 
is a complex number: A = |A|exp(i#). The order param- 
eter, as a function of coordinates, should be univalued. 
This implies that, as we travel along a closed contour, 
the phase of the order parameter can get an increment 
2mr where n is integer called winding number. The corre- 
sponding texture of the order parameter is called n-times 
quantized vortex. Let us now consider a singly-quantized 
vortex located at r = 0: 
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where <p is asimuthal angle in the x-y plane. The absolute 
value of the gap is assumed to be constant everywhere 
except regions of the size £ (correlation length) near the 
vortex centers. This assumtion constitutes the so-called 
London approximation. From Eqn.(Q) we can calculate 
induced charge density (per chain) using Eqn.(^J): 
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We see that it falls off as 1/r - it is highly non-local! 
However the screening makes the fall-off exponential. To 
see how this happens we make the following substitution 
in Eqn.(^|) in order to take vortex degree of freedom into 
account |]: (V0)k -► ik6 p k hason + {ik x z/k 2 )n k where 
n(r) — 27rYJni<5(r — r») is vorticity and Hi is the wind- 
ing number of a vortex at n . Minimizing the free energy 
functional with respect to gP hason we get an effective ac- 
tion: 
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Electric potential produced by vortices is obtained by the 
minimization of Eqn.(^) with respect to (p k and is given 
by: 
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C. Single Vortex 

Now n(r) = 27r<5(r) and the electrostatic potential pro- 
duced by a single vortex is given by the Fourier transform 
of Eqn.(|I^). It is hard to obtain an analytical expression 
in the general case, but one can get an idea of what this 
potential looks like from the consideration of the purely 
isotropic case: 7 = 1, Afc = A = const. With this simpli- 
fication the potential is given by: 
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Here J\ is first Bessel function. The effect of huge ac- 
tual anisotropy in A can be taken into account quali- 
tatively in the following way. First of all we may no- 
tice that there will be still two characteristic regions: 
r < e 1/2 A(6>) and r > e 1/2 A(#) where A(6>) is now angle- 
dependent crossover scale (which is, of course, extremely 
anisotropic). Then we may also notice that in the large- 
r region f(r) is still angle-independent: it is given by 
f(r) w A 2 v /r where A a „ is some "average" screening 
length. Important point here is that the total charge 
density in the large-r region falls off exponentially: one 
can check that by taking laplacian of ip(r). But the above 
mentioned extreme anisotropy comes to play in the small- 
r region - the core of a vortex is highly anisotropic and 
hard to analyze. 

Some of the vortex properties, however, allow exact 
description. From Eqn. (|l0|) it can be inferred that the 
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total charge of a vortex iz zero, but there is a non-zero 
electric dipole moment directed along y: 



d y = 2-KeeJX 2 = d 
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Here do = vp j2e. d y is exponentially increased as T — > 0. 
It should be noted that this expression was derived for a 
vortex in the bulk; it is expected that for vortices near 
boundary (the case which is relevant for vortices pro- 
duced by external fields - see below) the induced dipole 
moment is reduced. 



D. Critical Field 

Now we estimate the critical field at which an 
appearence of a single vortex becomes energetic ally 
favourable. In order to do that we substitute Eqn.dlC 
(with rik — 2n) back into Eqn. ([)]): 
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Here W is the size of the system in the y-direction. From 
Eqn.Q and Eqn.@ we obtain the critical field: 
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Comparing this expression with its superconducting ana- 
log one may call do an "electric dipole moment quantum". 
The dependence of the critical field on temperature is 
given through: A^ 2 = X t 2 fV / 2ttA /T exp(-A /T). So 
the critical field for vortex mixed state can be signifi- 
cantly lowered by increasing the screening length in the 
y-direction as T — > 0. As we do so, an interesting sit- 
uation can occur: D^ 1 can be made smaller than D\ 
signalling a temperature- driven transition between Type-I 
and Type-II-like ground states! A more intuitive expla- 
nation of this transition is obvious: one can, in principle, 
evaluate k x — X x /£ x and n y = X y /^ y . Usually one has 
k x -C 1 while IS exponentially temperature dependent. 
One can argue then that the type of the ground state is 
determined by the geometric mean of the two kappas: 
t\i — ■\ i J~^x^y • TlllS "mean" Ginzburg-Landau parameter 
k is also temperature-dependent. This raises a possibility 
of Type-I/Type-II temperature-driven transition as some 
critical value n c ~ 1 is passed by k in a temperature 
sweep. To our knowledge this is the first time that such 
type of transition is predicted. 



E. Screening Properties 

From the physical point of view the reason why the ap- 
pearence of vortices becomes favourable at high fields is 



that vortices can screen external electric fields. In order 
to show this we first introduce a "dense limit approxi- 
mation" H in which the everage distance between vor- 
tices is assumed to be much less than the characteristic 
length of the external electric field variation. Then one 
can take the vorticity n(r) to be a continuous function 
rather than the sum of delta- functions. The minimiza- 
tion of the effective energy Eqn.(^) with respect to 
becomes trivial: — d$ E% and leaves us with the ef- 
fective free energy which decribes the screening with the 
screening length A = Xtf- In fact this contradicts to 
the above made " dense limit approximation" , neverthe- 
less the conclusion about screening remains valid if we 
assume that the screening length is rather equal to the 
average inter- vortex distance: A = n~ x / 2 . But now the 
screening length becomes coordinate dependent leading 
to a non- exponential screening. Indeed, Poisson eqaution 
ip" — X~ 2 p = in this case reads ip" + d^ipip' — 0. Its 
solution is: 
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At distances ~ X qp , however, the exponential screening 
due to quasiparticles comes to play. So in the bulk the 
DW state is restored. Thus the electric field-induced 
transition to a vortex state is a surface effect. Therefore it 
is unlikely to be detected in any transport measurements 
p5| . On the other hand the best way to observe such 
a transition would be to measure directly the screening 
length. It would be expected that the screening length be 
~ X qp for field below critical whereas for fields above crit- 
ical it collapses to some small value. Such experiments 
are currently underway Jlq] . 



V. THERMALLY EXCITED VORTICES 

An alternative to the production of vortices in DW 
systems by applying transverse electric field at low tem- 
perature is their thermal activation in the temperature 
region near T c . We have already calculated the energy 
required to produce a singly-quantized vortex in the ab- 
sense of external electric field (Eqn.JIJ) - this energy is 
given by pi: 
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At low temperature this energy is big, but as we get closer 
to the transition point E a decreases rapidly because of 
decreasing rigidity K (it is proportional to the conden- 
sate density). Moreover, K is renormalized to smaller 
values in the vicinity of critical point by the fluctuations 
of the order parameter. The description of the dynam- 
ics of thermally excited vortices in external transverse 
electric fields is beyond the scope of this article and will 
be reported elsewhere; however a simple picture can be 
outlined here. The key point is that there is a gas of 
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thermally activated vortices at some finite temperature. 
We shall restrict ourselves to only two kinds of vortices, 
namely those with winding numbers rii = ±1 and dipole 
momenta d y = ±|dy| respectively. It is assumed that 
the excitation of vortices with higher winding numbers is 
suppressed be their bigger energies. Then applying trans- 
verse electric field one would be able to polarize vortex 
gas. The dielectric constant of the vortex gas in the case 
of weak fields is: 

Aimdl 

e = -jr* (18) 

Here n is the density of vortices, n also enters the ex- 
pression for the correlation length in the gas of vortex 
dipoles. Thus one would be able to extract n from ex- 
periments concerned with the response of DW materials 
to transverse electric fields. Such experiments are cur- 
rently underway jl6| . 

VI. CONCLUSIONS 

We have performed an extensive analysis of vortices 
in highly anisotropic quasi- ID DW systems. This work 
builds a more detailed and correct picture of the struc- 
ture and dynamics of vortices as compared to || and 
also predicts several unique properties of such vortices 
never encountered before with the case of their supercon- 
ducting counterparts. Among those are the possibility of 
a temperature-controlled phase transition between Type-I 
and Type-II-like phases; the non- exponential screening of 
strong electric fields by vortices at comparebly large scales 
as opposed to the incorrect conclusion about short range 
exponential screening drawn in Q ; temperature controlled 
"flux quantum" of a vortex - its dipole moment. At the 
end we would like to point out that the exploration of 
the dynamics of topological defects in DW systems is a 
relatively new and actively developing area of condensed 
matter physics with many surprises and, possibly, some 
connections to other branches of physics (see e.g. a recent 
paper jl7| and references therein). 
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